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Abstract 
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interactions are investigated. The approximate analytical expression for the ensemble-averaged 
spectral density for the Hamiltonian defined on the fully-connected graph is derived and analysed 
both for the electronic and vibrational problems which can be related to the contact process and 
to the problem of stochastic diffusion, respectively. It is demonstrated how to evaluate the extreme 
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I. INTRODUCTION 



The spectral properties of complex networks are of great current interest, both for practi- 
cal applications and from a fundamental point of view ,0,0,13. The 

physical phenomena occurring in the network can be described using the operators defined for 
the network. For example, the Hamiltonian describes electronic excitations in the network 
of atoms (nodes) characterized by energy levels communicating with each other by hopping 
integrals (links). The Laplacian operator for a set of atoms connected by elastic springs 
describes vibrational excitations [l3| or transport phenomena, e.g. stochastic diffusion of 
random- walk type 3,^^. The Liouville operator characterizes the spread of epidemics in 
;he network 16|, |l7| and the connectivity operator gives knowledge of the network topology 
n the matrix representation, these operators are characterized by matrices, the 
eigenspectrum of which gives rather complete information of e.g. dynamical properties of 
the network. The main difficulty in analytical evaluation of the network spectrum is in an 
inherent disorder incorporated into the network matrices. This can be topological disorder 
related to irregular Euclidean arrangements of the nodes (e.g. atoms in liquids and glasses) 
and/or disorder in connectivity (complex networks e.g. of scale-free or small- world type P|) 
or disorder in parameters associated with the nodes and links (e.g. mass and force-constant 

n 

disorder for the vibrational problem on a regular lattice llg or substitutional disorder in 
). 

The first task of the spectral analysis is to evaluate the spectral density of the relevant 
operator for a particular realization of disorder. For an ordered system, all the realizations 
are identical (zero disorder) and the spectrum of at least some operators can be easily found. 
For disordered systems, this is a highly non-trivial problem and can be solved implicitly only 
for some simple models of local perturbation in a reference not necessarily ordered system 
(see Ref. and references therein). The next task is to perform the averaging over dif- 
ferent realizations of disorder (configurational averaging). The configurationally averaged 
spectral density can be used then for comparison with the experimentally measured spec- 
trum (e.g. by inelastic neutron scattering for vibrational excitations [13]) or with other 
observable characteristics (thermodynamical values such as heat capacity or linear response 
functions, e.g. the dielectric response function). Such a comparison with experimental 
observables makes sense if the measurable quantity is self-averaging, i.e. if the difference 
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between the mean value of the observable for a particular realization of disorder in a macro- 
scopically large system and the cofigurationally averaged value of this observable tends to 
zero with increasing system size 2l|,|22]. Some of the observables such as thermodynamical 
characteristics of spin glasses away from the phase transition points and of normal glasses 
at not very low temperatures are self-averaging. However, some of the observables such as 



thermodynamical characteristics of spin glasses at criticality 23], low temperature electron 



conductance 



22] and dielectric response in disordered semiconductors with strong disorder 



2J| are sample- dependent and thus are not self-averaging. 



The spectral density of at least Hamiltonians and Laplacians exhibits decreasing fluctu- 
ations with increasing system size (even at the localization/delocalization threshold point) 
and thus is expected to be a self-averaging characteristic (see e.g. js^). This has been seen 
implicitly in numerous computer simulations of spectral properties of disordered systems (see 
e.g. 26| and references therein) and in experiment 27[. Usually, the ensemble averaging is 
tackled by means of mean- field theories ^ M with possible use of the replica method l3(| 
or introducing supersymmetry |2^|SH but in some cases the "exact" solutions are available. 
They are quite rare and the examples are the semicircular spectrum for the fully-connected 
graph (FCG) with random normally distributed node-node interactions 0, 3^ (see also 
Ref. j3J]) and Lloyd's model for a special type of the on-site disorder and any network 
topology 

The main aim of this paper is to present a model dealing with the disordered Hamiltonians 
defined on the FCG with a special type of the node-node interactions and a rather general 
type of the on-site characteristics. It is possible to find an implicit analytical expression for 
the spectral density for a particular realization of disorder and then to perform analytically 
the ensemble averaging of the spectrum of the Hamiltonian with precision up to 0{N~^) 
with being the number of nodes in the FCG and thus to demonstrate the self-averaging 
properties of the spectral density. The analytical results for the ensemble- averaged spectral 
density are available due to the existence of an exact solution for the matrix elements 
of the resolvent operator for a particular realization of disorder. The general solution is 
specified and analysed for several particular problems including the electronic and vibrational 
problems with multiplicative interactions defined on the FCG. The electronic problem is 
equivalent to the contact process in the dilute regime and the results can be used for the 
lower bound estimate of the critical point for the contact process which describes e.g. the 
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spread of epidemics through the network. The vibrational problem is equivalent to the 
stochastic transport problem and the results can be used for the investigation of dynamics 
of information packets propagating through a communication network. The main analytical 
results both for the electronic and vibrational problems are supported by direct numerical 
diagonalization of the Hamiltonian. 

The paper is organized in the following manner. The formulation of the problem is given 
in Sec. |n] Several simple examples are considered in Sec. Ill II followed by the general solution 
of the problem for polynomial interactions in Sec. IIVI The limitations of the approach are 
discussed in Sec. El The conclusions are made in Sec. ED and some derivations are presented 
in Appendices \M and IB] 



II. FORMULATION OF THE PROBLEM 



Let us consider a FCG containing nodes. Each node i is characterized by the parameter 
Ei (node bare energy) and the link between nodes i and j by parameter Vij (node-node 
interaction). Then we define an operator H ("Hamiltonian") on this FCG in the following 
manner, 

h = e(^^+^E^^^- + ^*0 ' (2-1) 

i \ jj^i / i,j 

where the self-interaction matrix element Vu is introduced for convenience (the Hamiltonian, 
in fact, does not depend on it). The tuning parameter 7 gives an opportunity to distinguish 
between two types of problems: (i) electronic-like for 7 = and (ii) vibrational-like when 



7 = 1 and all = (see also 



37[). For vibrational problem, the operator H is the Hessian 



operator and its elements obey the sum rule, Va = 'Ylj^i ^ijy which follows from the global 
translational invariance of the Hamiltonian |l3| |. 

Both the electronic and vibrational problems are usually defined on networks with Eu- 
clidean topology describing real materials. Below, we consider a FCG and thus the physical 
meaning of the Hamiltonian ()2.H) defined on the FCG should be specified. This can be done 
by introducing two mappings. 

First, the electronic Hamiltonian (7 = 0) is e quiv alent to the Liouville operator, C, 
describing the contact process in the dilute regime Indeed, the time evolution of the 
state vector, |-P(t)), for the contact process is governed by the master equation describing 
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the conserved probability flow dt\P{t)) = jC\P{t)), which can be rewritten in the dilute 
regime as (see Ref. |38| for more detail) 

dtP.it) = -riP,{t) + J2 Wj^Pjit) , (2.2) 

where Pi(t) is the probability of finding node i in an occupied state independent of the occu- 
pation of all the other nodes which can be in two states, occupied (infected) or unoccupied 
(susceptible), r, is the recovery rate for node i and Wji is the transmission (infection) rate 
between node j and i. The formal solution of Eq. (j2.2p is given by 

\m) = e^'im) = J2e^^'{e^m)\e^) , (2.3) 

j 

with Ej and |e-^) being the eigenvalues and eigenvectors of the Liouville operator, respectively, 
which coincides with the Hamiltonian ()2.H) for 7 = 0, = — Cj and Vij = —Wij. The long- 
time behaviour of the contact process in the dilute regime is defined by the maximum 
eigenvalue, ej^^nax, and if ej^max < then the epidemic goes to extinction and it invades if 
£j,max > 0. The approximate rate equation ()2.2|1 has been obtained by replacing the term 
{Pj{t) — Pjiif)) (where Pjiit) is the probability for both nodes i and j to be occupied 
independent of the state of all the other nodes) in the exact equation with Pj{t). Such an 
approximation enhances the transmission of the disease and thus the estimate of the critical 
point obtained from the solution of the equation, ej^max = 0, gives a reliable lower bound 
estimate of the critical parameter for the contact process meaning that if the disease does 
not spread in the dilute regime then it certainly does not spread in the system. This can be 
practically important for controlling epidemics in disordered systems where the estimate of 
the exact value of the critical parameter is a rather complicated task I39I. l4fl[ 41 1. 

The second mapping connects the vibrational problem to the problem of stochastic dif- 
fusion through a net. While the contact process describes a propagation of excitations 
(infected nodes) through the net with a not-conserved number of excitations (the num- 
ber of infected nodes changes with time) then the standard stochastic diffusion deals with 
propagation of the conserved number of excitations (diffusing particles) through the net by 
means of diffusional jumps (characterized by the rates Wij) between the nodes. The balance 
equation for stochastic diffusion coincides with Eq. ()2.2|) where rj is replaced by Xljyj 



mm 

which reflects the conservation of the number of particles. Therefore, under 



5 



the assumption of symmetric transition rates, Wij = Wji, stochastic diffusion through the 
network is described by the Hamiltonian ()2.1|) with 7 = 1, = and Vij = —Wij. For 
complex networks, such as the FCG, the diffusing particles can be associated e.g. with the 
information packets propagating through the communication network P]. The quantity of 
interest can be e.g. the return probability of the diffusing particle to the starting place, 
{Po{t)) = N-^Tt exp{Ht} (see e.g. |4J ). 

The aim of our analysis is to find the eigenspectrum of Hamiltonian ()2.H) defined on 
the FCG. In the site (node) basis, the Hamiltonian matrix is fully dense and, for the gen- 
eral case of arbitrary parameters and Vij, its diagonalization is not a trivial task and 
the solution is not currently known. However, for the classical case of a random matrix 
belonging to the Gaussian Orthogonal Ensemble, when the off-diagonal (diagonal) elements 
are independent and normally distributed (with variance doubled for diagonal elements), 
the configurationally averaged spectrum of semicircular shape can be evaluated analytically 
(with errors of 0{N~^)) (29l. l33l| . One of the key features of the matrices belonging to the 
Gaussian orthogonal ensemble is the statistical independence of the matrix elements. 

Below, we suggest another class of real symmetric matrices for which the spectral density 
can be found for a particular realization of disorder and then the ensemble averaging can 
be performed analytically. These are matrices with a particular (polynomial) type of the 
node-node interactions: 

Vi, = ipjacp^ , (2.4) 

where (pf is a n-dimensional row-vector, (pj = (1, (pi, (ff, . . . , 0"~^) and ck is a real symmetric 
n X n matrix of interaction coefficients, so that Vij{(j)i,(f)j) is a symmetric polynomial form 
of order 2n — 2 with respect to (pi and (pj, e.g. Vij = au + ai2{(pi + (pj) + a22(pi(pj for n = 2. 
The values of (pi are independent random variables characterized, in general, by different 
probability distribution functions, p^^{(pi). The order of interaction, n, is supposed to be 
much less than the number of nodes in the system, n N. The diagonal elements, Ei, in 
such matrices are also random independent variables distributed according to the probability 
distribution functions, PeX^i)- The probability distributions p^^i^cpi) and PeX^i) assumed 
to have all finite moments unless it is stated differently. 

We demonstrate below that the matrix elements of the resolvent (Green's function) op- 
erator, G, defined by the equation, (el — H)G = I, can be found exactly for the interactions 
given by Eq. ()2.4|) . This means, that the density of states (DOS), g{e), is available for a 



particular realization of disorder according to the following identity |45| . 

g(e) = -^lmTiGiE + iO) . (2.5) 

Due to the availability of the analytical expression for g{s), its configurational averaging, 
{g{e)), can also be undertaken analytically by means of the following integration (for N — > 

oo), 

{g{e)) = (^N-'J2^i^-^i)^^ J ■■■ I iV-'^^(^-^*)n/^<^«('^^)/^^-(^^)d'^^d^- (2.6) 
where the eigenvalues (eigenenergies) of the Hamiltonian. 

III. SIMPLE EXAMPLES 

Before considering the general case of the polynomial node- node interactions ()2.4j) . we 
present, first, three simple examples where the exact solution of the problem is available. 
These examples are for (i) the ideal FCG, (ii) the binary FCG and (iii) Lloyd's model 
defined on the FCG. In all cases, the biological (epidemiological) applications of the results 
are discussed. 

A. Ideal fully-connected graph 

A trivial case we need as a reference for further analysis is an ideal FCG, for which all 
the on-site energies are identical, PeX^i) = (^(ci — eo); and all the interactions are the same, 
P(l>i{.4'i) = ^{(pi ~ 4>o) (so that Vij = Vq), see Fig. HJ^a). The diagonal elements of the resolvent 
are, 

4°' = ^ + —^, (3.1) 

e - eo e - eo + NVq 

where eo = eo + Vq + 7(A^ — l)Vo, and thus the spectrum of the ideal FCG contains two 
delta- functions, one of them is (A^ — l)-degenerate, 

goie) = {gois)) = (l - {e - e~o) + ^H^ - + Wo) . (3.2) 

The spectrum of the Hamiltonian defined on the ideal FCG obviously obeys the "energy- 
conservation" principle, TrH = A^eo, meaning that the interactions do not change the total 
bare energy. If disorder is introduced in the ideal system it is quite natural to expect the 




(a) (b) 

FIG. 1: (Color online) Ideal (a) and binary (b) fully-connected graphs with N = 5. 

broadening of the (A^ — l)-degenerate level to the band and possibly the appearance of more 
levels split from the band. This is exactly what happens in the disordered FCG according 
to the analysis presented below. 

In terms of biological applications, Eq. ()3.2|1 can be used for estimating the value of 
the critical parameter rjc (this can be the ratio of the typical transmission and recovery 
rates) separating the absorbing (77 < ric) and active (77 > rjc) states of the system with 
respect to the spread of the contact process (epidemic). Indeed, the maximum eigenvalue 
of the spectrum coincides with the position of the non-degenerate 5-function and is equal 
to ^inax = {N — l)Wo — tq (bearing in mind that 7 = 0, Wq = — Vq and tq = —eg). The 
solution of the equation, ^max = 0, gives a standard mean-field estimate for the critical 
parameter 77* = (Wo/ro)* = {N — Obviously, if the transmission rate Wq is A^- 

independent then the critical value rj* approaches zero for large values of and the system 
is always in the active state 0]. However, if we assume that the transmission rate is inversely 
proportional to the number of nodes (in a migrating biological population the interaction 
time between the members of the population can be inversely proportional to the population 
size). Wo = wq/N, with wq being independent of A^, then the critical point exists and the 
estimate for the critical parameter is {wo/tq)* = 1. 
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B. Binary fully connected graph 



Another simple example of the network is a binary FCG (see Fig. db)), for which the 
spectrum for a particular realization of disorder is available and configurational averaging 
of the DOS can be performed exactly. The binary FCG consists of nodes of two types, A 
and B. The on-site energies and €3 and the node- node interactions Vaa, Vbb and Vab 
are defined by the types of the nodes. The node-node interactions, in general, are not of 
multiplicative form and can be described by Eq. ()2.4j) only if V^^ = VaaVbb- The only 
random parameter for the binary FCG is the number of nodes of a certain type, e.g. N^-, 
which is defined by the probability p (parameter of the model) for a node to be of type A. 
The values of A^^ or equivalently of concentration c = N^/N are distributed according to 
the binomial probability distribution, Pc(c) with the expectation value E[c\ = p and variance 
Var[c] = p{l — p)N~^ which is close to the variance of the normal distribution for N ^ 00. 

The spectral density for the electronic (7 = 0) Hamiltonian ()2.1|) defined on the binary 
FCG characterized by a particular value of concentration c contains four 5-functions, 

g(e) = (c - N-^)6{e - e^) + (1 - c - N-^)5ie - is) + " ' (3-3) 

j=l,2 

where = + Vaa and €3 = €b + Vbb- Only the two last 5-functions, 

g,(e) = N-' J2 6(6-6,) , (3.4) 

i=l,2 

depend on the random parameter c G [0, 1] and thus should be configurationally averaged. 
The values of 6, in Eq. ()3.4|) are the roots of the spectral determinant, D{6i) = 0, where 

V 6-€aJ\ e-es J [6-eA)[6-eB) 

The above expression for D{6) can be derived in a manner similar to the derivation for 
multiplicative interactions (see Appendix El) • 

For simplicity, we consider a symmetric binary FCG characterized by = = eo = — tq 
and Vaa = Vbb = Vq and also assume that the node-node interactions are negative and 
inversely proportional to A^, so that the interaction parameters, w = — A^Vq and A = 
—NVab, are positive and A^- independent. In this case, the c—6 map defined by the equation, 
D{6) = 0, with D{6) obeying Eq. ()3.5p . is given by the following bilinear form, 

(e + ro - w/2y + (A^ - w^){c - 1/2)^ = AV4 , (3.6) 

9 



where we have ignored the terms oc ^. Therefore, the positions of the S- functions in 
Eq. ()3.4p can be found as the roots of Eq. ()3.6|) for a particular value of c, 




ei,2 = -ro + -±\l--{X'-w^)i^c--j , (3.7) 

The form ()3.(j|l is hyperbolic for a weak interaction (A < w) between subgraphs A and 
B and elliptic for strong coupling (A > w). The configurational averaging in Eq. ()3.4|) is 
straightforward and 

/ / N\ N^^le + To - w/2\ r / N / M /o o\ 

^/(A^ - w^){X^/A -{e + ro- w/2f) 

where Ci,2 = 1/2 ± [(AV4 - (e + tq - w/2f)/{\^ - w'^)f^. The analysis of Eq. shows 
that, for any finite N , two (5-functions are broadened by disorder into two bands separated 
by a gap of width Wg = A for weak coupling and Wg = w for strong coupling (see Fig. IHl^a) 
and (b), respectively). The non-linearity of the map results in the singular behaviour of the 
ensemble- averaged DOS around e = — ro+w/2±A/2 where {gi{e)) oc {e+rQ—w/2^\/2)^^/'^. 

The ensemble-averaged spectrum for the binary FCG is bounded from the top by energy, 
£max, the knowledge of which is quite important from the applicational viewpoint (see below). 
The value of e^^x is easy to find from Eqs. ()3.6|) . ()3.8|) and it is 

£max = -ro + w for X<w, (3.9) 

and 

£max = -ro + ^(w + A) for \>w. (3.10) 

Eqs. ()3.9j) - ()3.10j) give the upper boundaries for the maximum eigenvalues for a particular 
realization of disorder, i.e. for a particular value of c. In the limit of large values of — oo, 
the binomial distribution Pc{c) is of the Gaussian form characterized by the negligible width 
and thus it approaches the ^-functional peak. Consequently, both energy bands collapse into 
two 5-functions located at e = e\^2 given by Eq. ()3.7p and ^max — Therefore, the critical 
value of parameter rj = w/tq can be found from the solution of equation, ei = 0. In the limit 
of weak coupling between the subgraphs, \ < w, the critical value rj* depends on c and lies in 
the interval, 1 < rj* < 2/{l + X/w). The value of t]* reaches the maximum, rj* = 2/(l + A/w), 
for the homogeneous FCG, when c = (FCG contains only nodes of type 5) or c = 1 (FCG 
contains only nodes of type A). The minimum, 77* = 1, is attained for equal concentrations 
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FIG. 2: (Color online) The c — e map (dashed line) for the symmetric binary FCG in the regime of 
weak (a) and strong (b) coupling. The ensemble-averaged DOS (scaled), {gi{e)), and the probability 
distribution function (scaled), Pc{c), are shown by the solid and dot-dashed lines, respectively. The 
values of the parameters are: (a) eg = 1, w = 1 and A = 0.1; (b) eQ = 1, w = 0.1 and A = 0.5. The 
probability p = 0.5 and = 100 for both regimes. 



of nodes A and B, i.e. when c = 0.5. In the regime of strong coupling between subgraphs, 
\ > w, the situation changes to the opposite one so that 2/(1 + X/w) < r]* < 1 and the 
lowest value of t]* corresponds to c = 0.5. 

Biologically, this means that mixed populations containing two species are most vul- 
nerable to epidemics for equal concentrations of species if these species strongly interact 
with each other (in terms of transferring disease) and if the species do not interact strongly 
enough then the mixture of species enhances the resistance of the population. 

The above estimates hold for the populations of organisms which are able to communicate 
(transfer a disease) to all other members of the population, i.e. for populations having 
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the communication network with the topology of the FCG. The other assumption is that 
the transmission rate between two individuals is inversely proportional to the number of 
individuals in the population. This is a plausible assumption for migrating individuals when 
the probability to establish a contact can be proportional to A^^^. 



C. Lloyd's model 

It is known that the configurational averaging of the spectral density can be performed 
analytically for the network of any topology, and thus for the FCG as well, if the diagonal 
elements of the Hamiltonian matrix are distributed according to the Cauchy distribution. 



6 

TC S"^ + (Ci - Co) 



where 5 is the width of the distribution, and all the relevant node-node interactions are not 
random, Vij = Vq (Lloyd's model Ssj]). 

The diagonal elements of the configurationally- averaged resolvent operator, {Gii{e)) can 
be expressed via the resolvent elements for the ideal FCG, G\'^\ with the argument shifted 
to the upper half of the complex plane, {Gii{e)) = Gfl\e + i6). Therefore the ensemble- 
averaged spectral density for Lloyd's model defined on the FCG has the following form: 

(gie)) = (l - ^) P^i^ - ^0 - l{N - l)Vo) + ^p.(£ - (7 - l)(iV - m) , (3.12) 

i.e. the DOS is obtained from that for the ideal FCG by broadening of both 5-functions in 
Eq. ()3.2j) into two Lorentzian peaks. Formally, the FCG with on-site energies distributed 
according to the Cauchy distribution is equivalent to the ideal FCG with nodes characterized 
by the complex on-site energies (see also 

In contrast to the binary FCG, the widths of both of the Lorentzian peaks in Eq. ()3.12|) 
do not depend on and for any value of 6 it is possible to find such a value of starting 
from which the positive eigenvalues appear in the spectrum. This means that Lloyd's model 
on the FCG is not resistant to the invasion of epidemics at least in the dilute regime. It is a 
consequence of the special form of the distribution of the recovery rates given by Eq. 1)3.1111 
with not existing high moments. 
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IV. POLYNOMIAL INTERACTIONS 



We have considered above several simple examples of the FCG for which the spectrum can 
be found exactly and its configurational averaging can be performed analytically. A natural 
question arises about the possibility of a similar analysis in the case of more general disorder. 
Below, we demonstrate that indeed in the case of the polynomial node-node interactions such 
an analysis is possible. 

A. General Solution 

In the case of the polynomial node-node interactions for the FCG, the matrix elements 
of the resolvent operator in the site basis can be found exactly and the spectral density for 
a particular realization of disorder is given by the following formula (see Appendix E} : 

. N Id 

g{e) = go{e) + 6g{e) = -J^^ie - e,) - — Im - [In (det D(e))] . (4.1) 

i 

Here the renormalized bare energies, are given by Eq. ()A2|1 and the spectral determinant, 
detD(£:), satisfies Eq. (|X7j) . 

Assuming for simplicity that all and (pi are identically distributed according to the 
probability distribution functions PeX^i) = P^i^i) P(t>i{'Pi) = P(f>{(pi), respectively, the 
spectral density g{e) can be configurationally averaged according to Eq. ()2.6|) . 

{g{e)) = (goie)) + {Sg{e)) 

OO 

/X d. 
P0(0i)Pe(e - Vu) - ^ ^ (arg (det T>{e))) , (4.2) 

— oo 

where 

N 

Vii = 7V'f "Pi + (1 - - iNif^I^ + (1 - -f)i/^fip^ , (4.3) 

j 

with 'i/)f = (fifct and ^ = J^^ipip^{(t>i) d(f)i so that Vu depends on the characteristics of 
node i only (see below for the justification of this approximation). The ensemble-averaged 
spectral density, {g{e)), has two contributions. The first contribution, {gQ{e)), is given by 
the convolution of two probability distributions, one of them, p^, is shifted along the energy 
axis due to the node-node interactions. If both distributions pe and are band-shaped (e.g. 
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normal, box or 5- functional distributions) then the function {go{e)) also has the shape of a 
band of typical width A defined via the widths of the distributions pe and P(f,. 

The other contribution to the spectral density comes from {6g{e)) and its magnitude 
is negligible in the main band region (due to the factor A^~^ in Eq. (j4.2p ) and finite out- 
side the band with {6g{6)) being in the form of several peaks (see below). The functional 
form of {Sg{e)) depends on the properties of the spectral determinant detD. It follows 
from Eq. ()A7|1 . that the spectral determinant is a random value which depends only on 
macroscopic sums, = aki = J2f 'Pii^ ~ ej)~\ i.e. det D = D{e] oq, ai, . . . , a2n-2)- Ac- 
cording to the central limit theorem, the values of are distributed around the mean 

1 /2 

value Ofc ^ NUki. in the peak region of width 5ak — {NYai[aki\) , i.e. the relative 
peak width of this distribution approaches zero in the thermodynamic limit (A^ oo), 
^cik/'o^ N~^/'^ 0. Bearing this in mind we can perform configurational averaging of the 
phase of the spectral determinant in Eq. ()4.2j) approximately (assuming, in fact, that the 
spectral density is a self-averaging quantity), 

(arg (det D(e;ao, a2n-2))) ^arg (det D(£:;aJ^,ar,-- -,0^;;^)) , (4.4) 

where 

oo 



ak{e) = N p,{ei)p^{(j)i)deid(j)i 

J J e - ei - Vii{(t)i} 

— oo 

oo oo 

= N [ (f)'^R{e-M<l)i))d<l), + iN I 0fj(£- i/,,(0,))d0, , (4.5) 



with 



oo 



R{z)=/ ^^dQ and liz) = -7ip,(z) . (4.6) 

Z 6j 



-oo 



Eq. (j4.4j) becomes exact for an infinite number of nodes. Note that the same arguments 
were used in the derivation of Eq. (j4.3|) for replacing J^f with resulting in Va being 
a function of (pi only but not other (pj for j i. 

Expressing the phase, arg (det D(e; oq, oT, . . . , a2„_2)), via real and imaginary parts of the 
spectral determinant, D = D + iD, and differentiating it with respect to energy according 
to Eq. ()4.2|1 we arrive at the following final expression for {g{e)), 

oo ^ ~ ' ^ / ~ 

{gie)) = {goie)) + {6g{E)) ^ [ p^(0,,)p.(e - Vu) d0, - ^ ^ ^~ > (4-7) 

-i D + D 
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where prime means differentiation with respect to e. This is the main result of the paper. 
Eq. ()4.7p allows the spectral density of the FCG with polynomial interactions to be calculated 
for rather general distributions of the bare energies, e^, and interaction characteristics, 0j. 
The functional form of the spectral determinant D{6) depends on a concrete formulation of 
the problem but it is irrelevant for the main band shape and can influence only the positions 
of discrete levels split from the main band. An alternative derivation of Eq. ()4.7|) by means 
of direct integration of the left-hand side of Eq. ()4.4|) (and thus demonstrating the self- 
averaging property for the spectral density) is presented in Appendix |Bl for multiplicative 
interactions in the electronic case. 

In the band region, where the contribution from (go^e)) is significant, both functions D{e) 
and D{e) are typically of the same order of magnitude, D D and D D ~D/A(ifn<^ 
A^), so that the contribution from {6g{6)) to the total ensemble-averaged spectral density is 
macroscopically small, {6g{e)) ~ (A^A)^^. This is not surprising and is a consequence of the 
particular form of the node-node interactions given by Eq. ()2.4|1 forcing the majority of the 
eigenvalues, Ei, of the Hamiltonian ()2.H1 (the roots of the spectral determinant, D{ei) = 0) 
to be bound between the consequent renormalized bare energies, i.e. Ei G (?j,?j+i). This 
property is similar to the well-known phenomenon of the spectral reconstruction caused by 



the interactions of one level 
(band) (see e.g. Refs. 



2G 



e.g. _ associated with the defect) with a continuum of levels 
36, 



Outside the main band region, where D ^ 0, the function {6g{e)) contributes in the 
form of the Gaussian peaks (see the explanation below), {6g{e)) = ^(-^[^«]) 
of width 6e^,i = ^/Vai\e^ and centred around the expectation value, E[e^,i] (the function 
A^(-E'[x], Var[x]) stands for the normal distribution of x with the expectation value E[x] 
and variance Var[a;]). The peak locations coincide with the roots of the real part of the 
spectral determinant, D(£'[e*j]) = 0, and the number of roots, m, cannot exceed the order 
of the spectral determinant, i.e. the order of interactions, m < n. Indeed, for energies far 
away from the band, £^/A ^ 1 where S = \e — e — Vjj\ (where Vjj = J Vjjp(f){(l)j)d(l)i and 
e = / ejPe(ej)dej), we can estimate the value of ~ N(j)'^/£, so that the spectral determinant 
is the n-th order polynomial of {N/S) with constant finite coefficients (under the assumption 
that the interaction coefficients aij do not depend on A^). It can have maximum n of A^- 
independent roots, so that \e^,i — e — Vjj\ oc A^. This means that the Gaussian peaks are 
mainly macroscopically separated from the band. 
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Therefore, the interactions of polynomial type ()2.4|) change the bare spectrum of the 
Hamiltonian ()2.1|) in the following manner: (i) the bare band is shifted and deformed and 
(ii) several isolated levels macroscopically separated from the band are formed. Such a 
picture will be supported below by detailed analysis of some simple cases of the low-order 
interactions for n = 2. 



B. Multiplicative interactions for the electronic problem 

We start the analysis with the simplest case of multiplicative (separable) node-node 
interactions, when the second-order interaction matrix contains only one non-zero element, 
a22^ so that 

V,, = (1 0,) I ° M I ^ I = a22<i>^4>j ■ (4.8) 



^22 








The spectral determinant for a particular realization of disorder is 



N 



D{e) = l + Y.7^^ (4.9) 

I 

with 

= + % -ei + ^a22N(l)4 + (1 - 7)a220i ■ (4.10) 

We consider below the electronic problem (7 = 0) mainly but the general results for arbitrary 
values of 7 will be presented when possible. 

The roots of the spectral determinant, bound between the renormalized bare 

energies, Ei G (ej_i,ej) (for 0:22 > 0) which is obvious from the functional form of D{e) 
given by Eq. (j4.9|) . Therefore the eigenvalues of the Hamiltonian are expected to be very 
close to the renormalized bare energies, Si = ei + 0{A/N), and the changes in the spectral 
density within the band region for ii due to interactions ()4.8p should be negligible (oc A^~^). 
This property is very similar to Rayleigh's theorem in the theory of vibrations in disordered 
systems (isf. 

The ensemble-averaged spectral density is given by Eq. ()4.7|1 with 



2 



D{e;a^) = l + Rea2iE) = 1 + N / a22(piR{e - Vu) d<Pi , (4.11) 
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and 

D{e;a^) = Im^ = / a220-/(£ - V^i) d0i , (4.12) 



where the functions R{z) and I{z) are defined by Eq. ()4.6p . The properties of the function 
{g{e)) have already been discussed above. In particular, in the main-band region the contri- 
bution from {6g{e)) is negligible (oc (A^A)~^) while outside the band it can exhibit Gaussian 
peaks. For multiplicative interactions, there is only one Gaussian peak outside the band, 
{g{e)) ^ {Sgi^)) — N{e^,a'^), with the peak position being the solution of the following 
equation D{e^) = 0, i.e. 

oo 

1 + iV y a22<plR{e, - Vu) d0, = . (4.13) 

— oo 

This equation can be solved approximately assuming that the level is split from the band 



far enough in comparison with the band width, i.e. 



- e - Vii 



/A » 1, 



~ e - Ara220' • (4.14) 

Indeed, we see from Eq. ()4.14|) that the distance between the isolated level and the main 
band is macroscopically large, oc (if the coefficient 0:22 does not depend on A^), which is 
consistent with the result for the ideal FCG (see Eq. ()3.2j) ) showing the similar structure of 
the spectrum with the (A^ — l)-degenerate level playing the role of the main band. 

In order to estimate the width of the Gaussian peak, we consider a particular j-th real- 
ization of disorder for which the random value of be estimated in a similar manner, 
Ej^, ~ A^~^ ti — ^"i- Fi'OKi this expression according to the central limit theorem, 
we conclude that the values of e^* are normally distributed with the expectation value given 
by Eq. ()4.14j) and with variance, 

^ ^ + Ara^,Var[0^ (4.15) 



This expression for the variance coincides with that given by Eq. ()B22|) more rigorously 
derived in Appendix El It follows from Eq. ()4.15|1 that the peak width, a, depends on A^ 
and it increases with N, a (x \/N, if the interaction coefficient 022 does not depend on A^. 
However, if a22 oc 1/A^, the peak width decreases with N, a oc and the peak collapses 

to a 5-function in the thermodynamic limit. 
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All the results presented above for the multiplicative interactions are supported in Ap- 
pendix El by an alternative derivation of Eqs. ()4.7|) . ()4.11|) - ()4.12|) for the ensemble- averaged 
spectral density using direct approximate integration of Eq. (j2.6|) . 

Numerical analysis confirms all the conclusions made above. Fig. |21 demonstrates both the 
main band (a) and the Gaussian peak (b) for the case of multiplicative interactions for the 
electronic problem (7 = 0) defined on the FCG. It is clearly seen how the interactions shift 
the Gaussian bare band (dot-dashed curve in Fig. El^a)) and change its shape. The isolated 
level is normally distributed (see Fig. Efb)) and macroscopically shifted down along the en- 
ergy axis. For all cases, the results of direct diagonalization are practically indistinguishable 
from those obtained in accord with analytical expressions. 

The evaluation of the integral in the expression for {gQ{e)) (see Eq. (j4.7|) ) has been 
performed numerically in the above example illustrated in Fig. El Such an integration 
becomes trivial if one of the probability distribution functions is a 5-function. For example, 
if the on-site energies are randomly distributed while the interactions are all the same, 
P(f>{(pi) = S{(pi — (po), the configurationally averaged spectral density coincides with the shifted 
distribution of the on-site energies, {go{e)) = Pe{£ — tt220o) (due to the linear map between 
ej and e, i.e. £ = Cj = Cj + a220o which follows from Eq. ()4.10|) ). On the other hand, the 
on-site energies can be all the same, Pe(ej) = 6{ei — eo), but the interactions are randomly 
distributed. In this case, due to the quadratic map between 0j and e, i.e. e = ii = 60 + ^220^ 
(see Eq. ()4.10|1 ). there are two contributions to the ensemble-averaged DOS from different 
branches of this quadratic map (see Fig. IH), 

1 



2Wa22{e - eo) 



+ P. 



(4.16) 



^22 / \ V ^22 

with e obeying the following inequality, a22(^— eo) > 0. The ensemble-averaged DOS exhibits 
a singular behaviour around the boundary of the spectrum at e ~ eo (see Fig. EI), 

^ , > ^ , (4.17) 

'a22[e - eo) 



if P0(O) is a finite non-zero value. From this example, we see that the bare spectral density 
can display quite drastic changes in shape depending on the type and probability distribution 
of the interaction parameters. 

In terms of epidemiological applications, the above analysis can be used for an estimate 
of the critical parameter r^c- Indeed, the low-bound estimate for ric can be easily obtained 
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FIG. 3: (Color online) The ensemble-averaged DOS (the main band, (9(e)) — {go{e)), in (a) 
and separate level, (5(e)) — {Sg^e)), outside the main band in (b)) for electron problem (7 = 0) 
defined on the FCG with multiplicative interactions. Both the on-site energies and interaction 
parameters (pi are normally distributed according to Pe(ei) = A^(0,0.09) (the dot-dashed line in 
Fig. H^a)) and p^{(j)i) = A^(l,0.09); 022 = 1- The dashed lines in (a) and (b) are plotted according 
to the first term in Eq. ()4.7() and Eqs. (|B21|) - HB22|) . respectively. The solid lines represent the data 
obtained by direct diagonalization of the Hamiltonian for = 1000 nodes and averaging over 10^ 
configurations. The DOS in Fig. Etb) is scaled by factor N. 

for the FCG with polynomial interactions by finding the position of the largest isolated 
eigenvalue, e^'^^, and solving the equation, e^^^^^r]*) = with e™'^'' being the largest root of 
the real part of the spectral determinant, D{e'^''^^) = 0. An important and general conclusion 
which follows from the analysis presented above is that the isolated normally distributed 
roots of the real part of the spectral determinant scale linearly with A^, i.e. e™'^^ oc if the 
interaction coefficients, aij, and the interaction parameters, (pi, and thus the transmission 
rates, Wij, do not depend on A^. This means that equation e'^'^^{ri*) = does not have 
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FIG. 4: (Color online) The (p — e parabolic map (dashed line) for the electron problem (7 = 0) 
defined on the FCG with multiplicative interactions (eo = 2, 022 = !)• The interaction parameters 
are normally distributed according to P(j,{4>i) = A^(l,0.16). The ensemble-averaged DOS (scaled), 
{g{e)) ~ {go{e)), and the probability distribution function (scaled), P(f){4>i), are shown by the solid 
(Eq. (|4.16|) ) and dot-dashed lines, respectively. The dots scattered around the solid line were 
obtained by direct diagonalization of the Hamiltonian (N = 3000) and by averaging over 300 
realizations of disorder. 

the solution independent of and the system does not exhibit the phase transition in 
the thermodynamic limit (it is always in the active state). In the opposite case, when 
the transmission rates are inversely proportional to A^, the maximum eigenvalue does not 
depend on and the transition exists at least for the Hamiltonian in the dilute regime. The 
concrete value of the critical parameter t]* depends on the particular type of the interactions 
and the probability distributions for recovery and transmission rates. 



C. Multiplicative interactions for the vibrational problem 

For vibrational problem defined on the FCG (7 = 1 and = 0) with multiplicative 
interactions, the ensemble-averaged DOS is given by Eq. ()4.7j) with PeX^i) = and 

Q = a22iV0^0 = , (4.18) 
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FIG. 5: (Color online) The double-log plot of {g{e)) ~ {go{e)) shown in Fig. 0] around the singu- 
larity. The solid line was obtained using Eq. ()4.17|) and the dots represent the results of direct 
diagonalization. The dashed line, oc e^^/^, is used as an eye guide. 

where /3 = a22N. Eq. ()4.18p reveals the linear map between 0j and e, i.e. e = ii = Pipicj), 
and thus for the main band, 




(4.19) 



If the coefficient 0:22 does not depend on then (3 oc N and the location of the main 
band {go{e)) also scales with A^. In the opposite case of P being independent of A^, the 
ensemble- averaged spectral density is A^-independent. Fig. IHl illustrates the high quality of 
the approximate expression for {go{e)) ~ {g{£)) in the band region (the solid line in Fig. IHl) 
by comparison with the results obtained by direct diagonalization of the Hamiltonian matrix 
(the dashed line in Fig. IH]). 

Similarly to the electronic problem, the contribution of {6g{6)) to the spectral density is 
negligible (oc 1/A^) in the main band region. A specific feature of the vibrational problem 
is that one of the peak-shaped contributions from {6g{e)) is always of the 5-functional form 
at exactly zero energy. For the multiplicative interactions, this is the only peak-shaped 
contribution, i.e. {6g{e)) = N~^6{e). This is a consequence of the global translational 
invariance of the Hamiltonian and also follows from the solution of the equation D{6^) = 0, 
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FIG. 6: (Color online) The ensemble-averaged DOS for the vibrational problem (7 = 1) defined 
on the FCG with multiplicative interactions. The random interaction parameters (pi are normally 
distributed according to p^{(t)i) = A^(l,0.09); (3 = 1. The dashed line is plotted according to 
Eq. 1)4. 19() and the data for the solid line are obtained by direct diagonalization of the Hamiltonian 
for N = 3000 nodes and averaging over 50 configurations. 

which reads 

00 

D(e,) = 1 + / ^'^j pM)'^<t^^ = • (4.20) 
J e^- I3(p(pi 
—00 

Obviously, the value of = is the solution of Eq. ()4.20|) and thus the peak is located at 
zero energy. This statement holds for an arbitrary realization of disorder, when the integrals 
in Eq. ()4.2()|1 should be replaced by finite sums and thus the peak associated with e^. is of 
5-functional form. 

Therefore, the vibrational problem defined on the FCG with multiplicative interactions 
has a relatively simple ensemble-averaged spectral density. It contains the main band which 
is obtained by rescaling of the probability distribution for interaction parameters and a 
5-functional peak at zero energy. 

The results obtained for the vibrational problem can be applied to the investigation of the 
transport properties of communication networks characterized by symmetric transmission 
rates Wij = —Vij > (i.e. /3 < 0) and Wy = Wji. For example, the above analysis for 
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separable node-node interactions is applicable for a network with communication channels 
charcterized by multiplicative functions 4^. In this case, the dynamical characteristics of 
a packet propagating through the network crucially depend on the functional form of the 
distribution, p^{(f)i), of the interaction parameters (pi. One of these characteristics is the 
return probability, (Pn(t)), of the packet to the starting point of its diffusion through the 



network (see e.g. 



oo 

{Po{t)) = N-'Tt exp{m} = j e''{go{e))de + ^ , (4.21) 

— oo 

where {go{e)) is given by Eq. (j4.19p . The spectrum of the Hamiltonian (j2.H) is negative 
except for one eigenvalue located exactly at zero which gives rise to the last term in Eq. 1)4.2111 
describing the random return of the packet to the origin. The time dependence of the return 
probability is dictated by the functional form of P(/,(0) in Eq. ()4.19p . If the distribution 
p^{4>) has the form of a band, (p G [^min, ^max], separated from zero by a gap, i.e. (prnin > 0, 
then the long-time behaviour of the first term in Eq. (14.21)1 for the return probability is 
exponential, {Po(t)) oc exp(— /900mint). On the other hand, if the distribution of starts 
from zero with p^{<f) ^ 0) oc 0" with a > — 1, the return probability decays with time 
according to the power law, {Po(t)) oc t^"^^ for large times. This property can be quite 
important for monitoring the localization properties of the information packets by means of 
choosing the appropriate distribution for the interaction parameters 0j. 



D. Polynomial interactions of higher orders 

To conclude this section about polynomial interactions defined on the FCG, we should 
emphasise that there are no principal difficulties in extending the above analysis to the cases 
of higher-order interactions with n > 2. The features of the ensemble- averaged spectral 
density will be the same as for the multiplicative interactions. The shape of the main band 
can vary significantly depending on the order of interaction and the probability distributions 
of the bare energies and interaction parameters. 

A typical example is shown in Fig. [3 for the electronic problem (7 = 0) defined on the 
FCG with polynomial interactions of order n = 6 characterized by a normal distribution of 
the interaction parameters and a (5-functional distribution of the bare energies. In this case, 
the convolution of the probability distributions in Eq. ()4.7)) is trivial and the band shape 
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is dictated by the non-linear (p — e map (see the inset in Fig. [7j) given by the polynomial 
of order 2n ~ 1 with the singularities in the DOS being due to the extremal points in the 
(p — e map. The ensemble- averaged DOS calculated according to Eq. (j4.7|) (the dashed line) 
is practically identical to the exact one (the solid line) obtained by direct diagonalization of 
the Hamiltonian. The isolated levels, ^j*, can be found if necessary by solving the equation 
for the real part of the spectral determinant, D{ei^) = 0. 




FIG. 7: (Color online) The exact (solid line) and approximate (the dashed line) ensemble-averaged 
DOS (over 10^ realizations) for the FCG (A^ = 3000) with node-node interactions defined by 
Eq. 1)2.4(1 in which n = 6. The elements of the Q-matrix are random values withdrawn from the 
uniform distribution defined in the interval, Oij S [—1,1]. The on-site energies and interaction 
parameters are distributed according to /Oe(e) = 5(e) and p0(</>) = A^(0,0.16), respectively. The 
non- linear (p — e map is shown in the inset. 
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V. LIMITATIONS 



In the previous section, we have demonstrated that the evaluation of the spectral density 
and its configurational averaging for the Hamiltonian defined on the FCG can be performed 
analytically in some special cases. Namely, this can be done for the binary FCG exactly and 
approximately for the FCG with a particular polynomial type of the node-node interactions. 
One of the restrictions for the polynomial interactions is that the order of interactions must 
be much less than the number of nodes in the FCG, n N. Generally speaking the solution 
given by Eq. ()4.7|) is valid for arbitrary value of n but it becomes "useless" for n > in 
the sense that the conclusions made above about the structure of {g{e)) do not necessarily 
hold for this case. For example, the maximum number of levels split from the main band 
due to polynomial interactions should be less or equal to the order of interactions, n, and 
therefore it can reach the value of if n > A^. This means that the main band described 
by {goi^)) in Eq. ()4.7p can disappear completely and a new band or set of levels arises due 
to the contribution from {6g{e)) (see an example below). The positions of levels split from 
the main band should be found by solving the n-th order polynomial for the real part of the 
spectral determinant. The complexity of this problem is not less than that of the original 
eigenproblem and thus solution ()4.7p becomes not very informative. 

The condition n <^ N can be broken for a very important type of the node-node interac- 
tions, Vij = V{\(j)i — depending on Euchdean distance, |0j — between nodes, where 
the (pi play the role of the node coordinates which can be random values (for simplicity, 
we analyse one-dimensional space). For example, let us consider a set of nodes randomly 
displaced from the sites in the ideal linear chain, so that (-E[0j] = ja; a is the mean distance 
between nearest sites and j is an integer). Assume also that the node- node interactions 
decay exponentially with the distance, 

V;, = Hexp{-[/3(0,-0,)]2} , (5.1) 

where j3 is the inverse interaction length and Vq is a constant. The interactions Vij given 
by Eq. (j5.ip can be expanded in a Taylor series generally containing an infinite number of 
terms and thus presented in the form of Eq. (j2.4j) with n ^ N. 

However, in the case of long-range interactions when the typical interaction length, 
is comparable to the system size, L, i.e. (3^^ ^ L = Na, the Taylor series for Vij contains 
only a finite number of terms, n N, and the interactions are of the polynomial type. 
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Consequently, the ensemble- averaged spectral density should be well approximated in the 
main band region by the function {go{e)) evaluated according to Eq. ()4.7|) which gives for 
the electronic problem, 

(goie)) = p,{e - Vo) , (5.2) 

because ii = €{ + Va = ei + Vq. It follows from Eq. ()5.2j) that {go{e)) does not depend on the 
distribution of (pi at all. Indeed, we have found numerically that in the case of long-range 
interactions the ensemble-averaged spectral density follows the theoretical prediction ()5.2|) 
(the black curve in Fig.[7|obtained by direct diagonalization for j3Na = 1 is indistinguishable 
from that obtained using Eq. (|5.2p ). In this regime, all the nodes interact with each other 
at approximately the same strength, Vij ~ Vq, and the system is equivalent to the FCG with 
on-site energy disorder only. 

When the typical interaction length decreases, more terms in the Taylor series for Vij 
should be kept for the accurate representation of function ()5.1|) leading to an increase in the 
order of the polynomial interactions. An increase in the value of n gives rise to more and 
more levels split down off the main band (see the red solid curve in Fig. IHlfor (3Na = 10^). 
These separate levels eventually form a broad band for medium-range interactions (see the 
dark green dashed curve in Fig. |H1 for f3Na = 5 x 10^) which transforms to a relatively 
narrow band of width, A ^ 4Vo exp(— /5^a^), having a well-recognizable shape of a band 
for the spectrum of an ideal linear chain with nearest-neighbour interactions broadened by 
diagonal and off-diagonal disorder (see the blue double-dot-dashed and light green solid 
curves in Fig. IHlfor jSNa = 10^ and jSNa = 1.5 x 10^, respectively). The original main band 
centred around e ~ e + Vq for long-range interactions eventually disappears with increasing 
value of p. 

Therefore, we can conclude that the theory for polynomial interactions presented in 
SecHVlcan be applied to systems with Euclidean long-range interactions but fails to describe 
the short-range interactions. 

The other limitation of the above analysis concerns its applicability to the disordered 
complex networks of the FCG topology only. Of course, the topology of real complex 
networks is much more complicated (e.g. scale- free or small- world topologies j^) and the 
FCG can be considered as a first approximation for the networks with high node-node 
connectivity. 



26 



3- 



V 



2 - 



1 - 













- M-i 












■ P-i 

P=1.5 





-2 




FIG. 8: (Color online) The exact ensemble-averaged DOS (over 3 x 10^ configurations) for a linear 
chain on = 10'^ nodes interacting with each other according to Eq. (jS-lf) (Vq = 1 and a = 1) for 
different values of the inverse interaction length /? (as marked). The on-site energies and interaction 
parameters are normally distributed according to Pe{e) = A^(0, 0.01) and p^{cj)) = A^(l,0.01), 
respectively. 

VI. CONCLUSIONS 



To conclude, we have analysed the spectral properties of the Hamiltonian both for elec- 
tronic and vibrational problems defined on the fully-connected graphs with a special type 
(polynomial) of interactions. 

Our main finding is the analytical formula (see Eq. I4.7|l for the ensemble-averaged spectral 
density. The ensemble-averaged spectral density has two contributions with clear physical 
interpretations: (i) the first contribution describes the main spectral band and (ii) the 
second one is related to the set of discrete levels separated from the band. The main band 
originates from the bare spectral band shifted and deformed by means of a convolution 
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of two probability distributions of bare energies and interaction parameters. The discrete 
levels are split from the bare spectral band due to interactions and the number of such levels 
depends on the order of the interactions. 

The approximate analytical configurational averaging is possible due to the availability 
of the exact analytical solution for the resolvent matrix elements for a particular realization 
of disorder (see Appendix^. Technically, configurational averaging is done with the use of 
the central limit theorem by replacing the function of fluctuating macroscopic sums by the 
function of the means of these sums (see Eg. 14. 4|) . This step can be justified and illustrated 
by a more rigorous approach for a particular case (see Appendix El) • AH the final results are 
convincingly supported by numerics. 

Both the electronic and vibrational problems discussed in the paper can be mapped to 
the contact process in the dilute regime and the stochastic diffusion problem, respectively. 
The results for the electronic case can thus be used for obtaining the lower bound estimate 
of the critical parameter for the contact process while the results for the vibrational case 
allow the return probability for an information packet to be evaluated for the communication 
networks. 

APPENDIX A: EXACT SOLUTION FOR THE RESOLVENT MATRIX ELE- 
MENTS 

The resolvent matrix elements for the Hamiltonian given by Eq. ()2.ip can be found exactly 
for the node- node interactions of polynomial type ()2.4j) . In order to demonstrate this let us 
recast the equation for the resolvent, {el — H)G = i, as 




where the renormalized bare energy, ej, is 

N 

e, = e, + Yl + (1 - = ^i + Vii , (A2) 

j 

with ipf = (pjcx. The value of Vu = •yil^f J2f + (1 ~ l)'^J'Pi = ^i + Vu generally depends 
on all (f)j but Va = Vij(0j) for the electronic problem (7 = 0). 

The node-separable type of polynomial interactions in Eq. ()2.4j) {i — j interaction is 
proportional to the matrix product of separate node characteristics) allows factorization to 
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be performed there by introducing xj = X^kVkGkj, so that 

G., = -ZT - J^T^-^J • (A3) 
Eq. ()A1|) can be muhiphed by ^pi and summed over i and thus transformed to equation 
which can be solved for xj, 



Substitution of xj into Eq. ()A3|1 resuhs in the required result for the resolvent matrix ele- 
ments, 

G,, = ^-^{l + y^\ (A5) 
The DOS is expressed via TrG, which is given by the following expression, 

T^6 = EF^-Ei^D-'^, (A6) 

and 

U.{,,t^y (AT) 

The eigenvalues, e^, of the Hamiltonian coincide with the poles of the resolvent which are 
the roots of the spectral determinant, D{e) = detD(£:)|, i.e. D{ei) = 0. This follows from 
the form of Eq. ()A6|) in which the contributions to the denominator from the first sum, 
oc {e — ii), are cancelled by the similar terms in Eq. ()A6j) . 

It is easy to show that Eq. ()A6|) can be recast in an elegant form via the derivative of the 
spectral determinant. 



N 

e — li ' de 

Indeed, rewriting the second term from Eq. ()A6|) in the following form 



^ 1 dln[det(D)l , , , 

TrG = > + \^ . A8 



f 4Ld-'^ = Idet(D)]-' f "-^-cj^ , (A9) 

I k,m 



where C stands for the matrix of cofactors for matrix D, and comparing it with the derivative 
of det(D), 

m,k m,k 

we arrive at Eq. ()A8|) . Finally, taking the imaginary part of Eq. ()A8|) and using then Eq. ()2.5p 
leads to Eq. (gH). 
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APPENDIX B: CONFIGURATIONAL AVERAGING BY DIRECT INTEGRA- 
TION 



In this Appendix, we give an alternative derivation of the expression for the ensemble- 
averaged spectral density in the case of the multiplicative node-node interaction, Vij = 
a22<Pi(t>ji for the random Hamiltonian defined on the FCG. The derivation is similar in some 
aspects to that given in Ref. [2^ for the mean density of eigenvalues of a 2D integrable 
billiard. 

The starting point for the derivation is Eq. ()A6j) recasted for the multiplicative interac- 
tion (|4.8|) in the following form: 



N ^ N ,9 / N 



0- 



l .y^ a220- f _iMi+Ef 

e -ii ^ (e- h 



where 



ej = ej + a220i , (B2) 

for the electronic problem (7 = 0) analysed below for concreteness (the analysis can be 
easily extended to the vibrational problem (7 = 1)). Using definition (|2.5|) and Eq. (|2.6p we 
obtain the expression for {g{e)) = {go{s)) + {6g{e)), where {goi^)) coincides with the first 
integral term in Eq. ()4.2|) with Va = a220f and {6g{e)) is given by 

0~ N 12 2 

me)) = I e'^dk J-. -J ^^^^ J]e''=Wp^(0,)p,(e,)d0,de, , (B3) 

or by the equivalent expression (assuming for definiteness that 0:22 > 0), 

0- 

{5g[e)) = -Re [ e"" [F{k, e)f~^ Q{k, e)dk , (B4) 



TT 



with 



00 



F{k,e)= I p^{(f)i)f{(f)i,k,e)d(f)i= I p^{4)i) I e'' -'^^^ p,{ei)deid(j)i , (B5) 



and 



Q{k,e)= I p^{(f)i)q{(f)i,k,e)d(j)i = I p^{(j)i) I z^e p,{e,)de4(j), , (B6) 
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where the function k,e) is related to f{(pi, k,e) by the following equation, 



a22(pj ok^ 



The next step is in the evaluation of the integral, 

f{(f)i,k,e)= / e"'^p,(ei)dei , (B8) 



{x = —ka22<f)'i > 0) with the integrand exhibiting an essential singular point at ej = e + iO — 
a220i = z. This can be done by expanding the exponential function in a Taylor series and 
integrating each term. 



/= / Pe(ei) + ^^7 

z — e,- 2 [z — e 



de,; 



1 - ix {R{z) + il{z)) + ^ {R'{z) + U\z)) + O(xVA') , (B9) 

where the functions R{z) and I{z) are defined by Eq. ()4.(jj) . In order to evaluate the third 
term in Eq. ()B9j) we integrated once by parts and used the following identity, 

CO oo 

= / M^l^d.. = / ^d.. . (BIO) 

J [z-eiY J z-et 

— oo — oo 

In expansion (|B9|) , we keep only the terms up to the second order in /?/ A including because 
this is enough for obtaining the leading term for k^e) according to Eq. ()B7|) 

k, e) = -a220- [R'{e - a220- ) + il'{e - a220' )] • (BH) 

The higher order terms in /5/A are small both for f{(f)i,k,e) and q{(f)i,k,e) because typ- 
ical values of \k\ significantly contributing into integral ()B4|) in the band region are 
\k\ ~ A/{a22<P'iN) (see below) and thus (3/ A ~ ^ 1. Note, that the same result 
for f{(j)i,k,e) can be obtained by a different method based on the shift of the essential 
singularity to infinity as was suggested in Ref. [2^. 

Substitution of Eqs. (fRTn|l - (fRTT^ into Eqs. ^^-^^ gives, 

F{k,e) = l + iA;[/?i(£)+i/i(£)] + y [i?^(e) + i/^(£)] (B12) 
Q{k,e) = -R[{e)-U[{e) , (B13) 
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where 



-oo 
oo 



(B14) 



(B15) 



Im{k,e) = J p^{(j)i) [a22(f)t]"' I{e - a22(f)'^i)d(f)i . 

— oo 

Using the above expressions for F{k,e) and Q{k,e) we can rewrite Eq. ()B4|) as 

(^g^s)) = -^Rel[R[{e) + U[iE)]Je^^''^dk} , (B16) 



with 



ip{k) =ik + {N -1) In 



l + k{iRi-Ii) + —iR'^ + il!, 



= iA;(l + NR^) - kNh + ^ {RI + R'2 - /i) + {I2 + ^RJi) + 0{k^) . (B17) 
There are two energy regions: (i) inside the band where Rmi^) ~ Imi^) ~ (^22^"^/^ ^^d 



fii) outside the band where Rmie) 



a 



— el with le* — el ^ A and /mfe) either 



approaches zero (e.g. exponentially for the normal distribution p^) or identically equals zero 
for the box distribution. In these regions, integral ()B16|) has different contributions to the 
total ensemble-averaged spectral density. Inside the band, we can ignore the terms oc fc^ 
in expression ()B17|) for ip{k), so that e'^^''^ exponentially decays for |A;| ^ cxd on the typical 
scale k ~ l/(A^/i) ~ A/(a220^iV), and 

1 (1 + NRi)NI[ - NR[NIi 



(He)) 



ttN (1 + NRiY + {Nil 



(B18) 



which exactly coincides with the second term in Eq. ()4.7|) bearing in mind that D = 1 + NRi 
and D = Nil. 

Outside the band, ^ and thus 0, so that the real linear term in k in Eq. ()B17|) 

becomes negligible and the next terms in the expansion must be kept. 



k'^N 

ij{k) ~ iA;(l + NRi) + ^[Rl + R'^) 



where 



a220i 



e - a22<Pi 



;)pe(ei)d0idei = -x^ 



(B19) 



(B20) 
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and Ri = x,Xi = a220?/(£^ — «220i — q)- Note that Eq. (IB19|1 contains an extra term, oc Rf, 
in comparison with a similar expression given in Ref. j2o| which is due to a more accurate 
expansion of ^/^ in /c. Straightforward evaluation of the Gaussian integral in Eq. ()B18|) leads 

to 

{Sg{6)) = ^ e 2ivva.-K] = — e 5^ , (B21) 



v/27rA^Var[x,] ^2 



2 



with Var[xi] = —Rl — R'2 and 



= - ^ (Varh] + Var[iVa220?]) • (B22) 



The spectral density given by Eq. ()B2Hl represents a Gaussian peak of width a centred at e = 
The location of the peak ~ e — 022^^ is the solution of the equation, 1 + NRi{e^) = 
identical to Eq. ()4.13|) and thus the expression for e^, coincides with that given by Eq. ()4.14|) . 

All the derivations presented in this Appendix were undertaken under the assumption 
that the coefficient 0:22 does not depend on A^. However, in the case when 022 oc all the 
results obtained for the main band still hold but Eq. ()4.14|1 for the position of the separate 
levels is no longer correct in general and the equation, 1 + NRi{e^) = 0, should be solved 
without using the assumption that the level is well separated from the main band. 
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